Mechanics of bundled semiflexible polymer networks 
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While actin bundles are used by living cells for structural fortification, the microscopic origin of the 
elasticity of bundled networks is not understood. Here, we show that above a critical concentration 
of the actin binding protein fascin, a solution of actin filaments organizes into a pure network of 
bundles. While the elasticity of weakly crosslinked networks is dominated by the affine deformation 
of tubes, the network of bundles can be fully understood in terms of non-affine bending undulations. 

PACS numbers: 87.15.-v, 87.15. La 



The mechanical properties and dynamic organization 
of the cytoskeleton determine the morphology and me- 
chanical response of eukaryotic cells. To ensure adapt- 
ability of both organization and mechanics cells exploit 
the dynamic interplay between semi-flexible polymers 
such as microtubules and actin filaments using a mul- 
titude of associated binding proteins. In particular, the 
local elastic properties are regulated by the activation of 
auxiliary proteins which e.g. cross-link and/or bundle the 
filamentous networks into complex scaffolds. Given the 
importance of the actin cytoskeleton for force generation 
and transduction there is much interest in understanding 
the mechanical properties of different network structures 
and the physical origin of the transitions between them. 
This is best studied in in vitro model systems [l[ . In the 
absence of cross-links actin solutions are successfully de- 
scribed by the spatial confinement of thermal bending un- 
dulations upon affine tube deformation 0). Cross- linked 
semiflexible polymer networks, on the other hand, are 
in general dominated by an interplay between polymer 
stretching and bending modes, the precise form of which, 
as well as the degree of non-affinity, strongly depends on 
the network microstructure 0] . So far the mechanical re- 
sponse of highly cross-linked actin networks, also in the 
presence of bundles and composite phases, has mainly 
been described assumingpurcly affine cntropic stretch- 
ing deformations 0,0, 0,13 • However, an applied tension 
can stretch a thermally undulating polymer only as far 
as there is excess contour length available. As the maxi- 
mal amount of stored length is inversely proportional to 
the persistence length, entropic stretching is suppressed 
in networks of stiff polymer bundles, where the persis- 
tence length grows with bundle size [H, Moreover, 
the highly non-linear nature of the force-extension rela- 
tion of semi-flexible polymers implies that linear elastic- 
ity is applicable as long as only a fraction of the total ex- 
cess length is pulled out. As an alternative the recently 
introduced concept of the "floppy modes" may be bet- 
ter suited to describe the polymer elasticity in situations 
where entropic effects arc suppressed These floppy 



modes constitute bending excitations which, unlike the 
affine stretching deformations, retain a highly non-affine 
character. 

In this Letter we show that above a critical concentra- 
tion of the actin binding protein (ABP) fascin, a solution 
of actin filaments organizes into a homogeneous network 
whose building blocks are bundles only. At low cross- 
linker concentration, the network response is dominated 
by the affine deformation of reptation tubes and the en- 
suing changes in confinement free energy The ob- 



served mechanical and structural transition between both 
phases can be described by a simple relation between the 
ABP concentration and the entanglement length. It is 
proposed to rationalize the scaling of the elastic modulus 
in the bundled regime in terms of the floppy mode pic- 
ture. A model based on affine stretching deformations 
only fits the data if additional assumptions about the 
bundle structure are made. 

G-actin is obtained from rabbit skeletal muscle and 
stored in lyophilized form at -21 °C [l^ . For measure- 
ments the lyophilized actin is dissolved in deionized wa- 
ter and dialyzed against G-Buffer (2 mM Tris, 0.2 mM 
adenosine triphosphate (ATP), 0.2 mM CaCl 2 , 0.2 mM 
dithiothreitol (DTT) and 0.005 % NaN 3 , pH 8) at 4 °C. 
The G-actin solutions are kept at 4 °C and used within 
seven days of preparation. The average length of the 
actin filaments is controlled to 21 um using gelsolin which 
is prepared from bovine plasma serum following [l3| . Re- 
combinant human fascin (55 kD) was prepare d by a mod- 
ification of the method of [l4[ as described by [l~5| . In the 
experiments the molar ratio R between fascin and actin, 
R = Cf /c a , is varied over almost three decades. 

To resolve the structure and mechanical properties of 
actin/fascin-networks actin is polymerized in F-buffcr 
(2 mM Tris, 2 mM MgCl 2 , 0.2 mM CaCl 2 , 0.2 mM DTT, 
100 mM KC1 and 0.5 mM ATP, pH 7.5). For fluores- 
cence microscopy filaments are stabilized with tetram- 
ethyl rhodamine iso-thiocyanate (TRITC) phalloidin; ei- 
ther labeled reporter filaments (1 per 400) or continuous 
labeling is used at distinct amounts of fascin. To avoid 
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photoblcaching 0.6 pM Glucoseoxidase, 0.03 pM Cata- 
lase and 0.01 M Glucose are added. The samples for 
transmission electron microscopy (Philips EM 400T) are 
adsorbed to glow-discharged carbon-coated formvar films 
on copper grids. The samples are washed in a drop of dis- 
tilled water and negatively stained with 0.8 % uranyl ac- 
etate; excess liquid is drained with filter paper. The vis- 
coelastic response of actin/fascin- networks is determined 
by measuring the frequency-dependent viscoelastic mod- 
uli G'(u>) and G"(uj) with a stress-controlled rheometer 
(Physica MCR 301, Anton Paar, Graz, Austria) within a 
frequency range of three decades. Approximately 520 pi 
sample volume are loaded within 1 min into the rheome- 
ter using a 50 mm plate-plate geometry with 160 pm plate 
separation. To ensure linear response small torques are 
applied. Actin polymerization is carried out in situ, mea- 
surements are taken 60 min after the polymerization was 
initiated. 

Fluorescence images show that in the presence of high 
concentrations of fascin, actin filaments organize into a 
network of bundles (Fig. QJi) while below a critical value 
R* ps 0.01 no bundles can be observed. Both fluorescence 
and transmission electron microscopy do not show any 
sign of composite phases or microdomains as observed 
in the presence of other ABPs @, 0, 17 1. Moreover, 
the existence of a purely bu ndled phase is demonstrated 
by a cosedimentation assay [Supplement | . The bundles 
formed are very long (> 100 pm) and straight which 
is consistent with the measured bending rigidity K [1]. 
TEM micrographs reveal that above R* the actin/fascin 
bundle thickness D and therefore the number of actin fil- 
aments per bundle, N, increases weakly with R (Fig.[Tp). 
The bundle thicknesses are extracted from the TEM- 
micrographs by fitting a Gaussian to the intensity pro- 
files, obtaining a scaling of D ~ N 1 / 2 ~ R x with 
x = 0.27. Concomitant with the structural changes 




FIG. 1: a) Fluorescence micrograph of an actin/fascin net- 
work (0.1 mg/ml actin): for high fascin concentrations a 
purely bundled network is formed (scale bar is 10 pi), 
b) From TEM pictures (inset, scale bar is 0.2 pm) a scaling 
relation for the average bundle diameter D is obtained. 

the viscoelastic properties of the network alter: with in- 
creasing R both the storage modulus G'(uj) and the loss 
modulus G"(lu) increase over the whole frequency range 
probed. The storage modulus G'{uj) exhibits a plateau 
at low frequencies, while the loss modulus G"(uS) reveals 
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FIG. 2: Plateau modulus Go as a function of the molar ratio R 
of fascin with respect to actin for two different concentrations 
of actin: 0.4 mg/ml (circles) and 0.2 mg/ml (squares). The 
dependence of Go on c a is obtained by scaling the fits for the 
0.2 mg/ml actin data upon the 0.4 mg/ml data points. The 
dashed line shows the boundary separating the two scaling 
regimes. The original frequency spectra for 0.4 mg/ml actin 
at different cross-linker concentrations (R = 0, 0.001, 0.002, 
0.005, 0.01, 0.02, 0.05, 0.1, 0.2, 0.5) are depicted in the inset. 



a well-defined minimum which shifts to higher frequen- 
cies with increasing R. The plateau modulus G'(10 mHz) 
plotted against R shows two distinct regimes in the elas- 
tic response. At low R, Gq is only slightly dependent 
on R, G ~ R a l ± °'\ while above a critical value R**, 
G increases with G ~ R 15 ± 02 (Fig.©. This expo- 
nent fits the data for both actin concentrations probed 
(c a = 0.2 mg/ml and c a = 0.4 mg/ml). The transition 
point R** agrees well with the structural transition at 
R* observed in microscopy. Below R* = R** the plateau 
modulus scales with the actin concentration as Go ~ c a 3 
suggesting that entanglements dominate the elastic re- 
sponse [2| . Above R* a different scaling regime occurs 
with G ~ c 2 a A . 

With the observed scaling behavior Gq(R, c a ) the 
plateau modulus is parameterized in both regimes, be- 
fore and after the structural transition. At the cross-over 
concentration R = R* these two parameterizations have 
to be equal. This uniquely determines the scaling of R* 
with the actin concentration, R* ~ c~ 79 , which results 
in the constraint Cf ■ c^ 21 ~ 1. This can be approx- 

1 /2 

imated to c/ ■ l e ~ 1 using the entanglement length 

— 2/5 

l e ~ c a ■ This surprisingly simple criterion for the 
bundling transition defies an obvious explanation and a 
detailed theoretical model is still lacking. It would need 
to account for the subtle interplay between confinement 
free energy of polymers in both the bundle and the net- 
work as well as the binding enthalpy of the cross-linking 
proteins. 
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The mechanical properties inside the bundled regime 
may, on the other hand, be understood in terms of the 
non-affine floppy mode model [T3| , where network elastic- 
ity is attributed to bending modes of wavelength compa- 
rable to the distance between cross-links, l c , and with 
stiffness fej_ ~ K>/lc- I n this picture typical deforma- 
tions of the network do not follow the macroscopic strain 
affmely but scale as <5 na ~ jLb, where Lb is a constant 
length over which an individual bundle within the net- 
work can be assumed to be straight. From our fluores- 
cence and TEM pictures, we would expect this length to 
be comparable to the bundle length. As a consequence 
the linear clastic modulus reads 



Go ~ vk±_5. 
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with the polymer density v ~ 1/£ 2 Z C . This model can 
be tested by relating the structural parameters of the 
network, mesh size £ and l c , and the bending elasticity k 
of the bundle segment to the concentration of actin and 
fascin monomers (c a ,c/). 

The structural information obtained by TEM and flu- 
orescence microscopy justifies the assumption that the 
bundles form an isotropic network similar to an entan- 
gled structure of single filaments. With increasing R, 
filaments and smaller bundles reorganize to form larger 
bundles that are spaced further apart. The mesh size £ 
of this self-similar structure therefore depends on R as 
£ ~ Co^V 1 / 2 , where £ ~ Ca ^ 2 is the mesh size of the 
filamentous network. Cross-linking will typically occur 
on the scale of the entanglement length l e , which plays 
the role of a distance between bundle-bundle intersec- 
tions (entanglement points). Since on average only a 
fraction of those will be occupied we can assume that 
distances between cross-links along the same bundle are 
given by l c ~ R~ v l e @, Q. Doubling the cross- linker 
concentration R should halve the distance between them, 
suggesting an exponent y m 1. 

For a description of the elastic properties of the bun- 
dles it is necessary to realize that fascin only leads to 
loosely coupled bundles, where bending is dominated by 
the shear stiffness of the cross-linking proteins [H, Q . The 
key quantity in this context is the bundle coupling pa- 
rameter a(l c ) = (l c /b) 2 , where the length-scale b ~ 5 1 / 2 
encodes the properties of the ABP's inside the bundle, 
in particular via the average distance S between cross- 
links. In general, S will depend on the concentrations c/ 
and c a , however, the precise relationship is not known. 
From fluorescence images the mesh size of the bundled 
network at R = 0.5 can be approximated which allows 
to calculate l e and thus l c . From this one can estimate 
the coupling parameter to be a > 1 for the whole bun- 
dle regime, implying that the effective bundle bending 
stiffness K acquires a wavelength dependence , leading 
to k(X) ~ NKfCt(X), where A is the wavelength of the 
deformation mode. This stands in marked contrast to 



what is known for single filaments or scruin-bound bun- 
dles where a fully coupled bending regime, k ~ N 2 Kf, 
has been assumed [18| . The wavelength dependence of 
the bundle stiffness has far reaching consequences for the 
static as well as dynamic properties of semi-flexible poly- 
mer bundles In particular, it implies that the en- 
tanglement length l e has to be reevaluated. As it derives 
from the suppression of long wavelength fluctuations by 
confining a polymer into a tube it is highly sensi- 
tive to a wavelength dependent k(A). This results in 
~ (Nip)!; 4 /b 2 , which is different from the usual ex- 
pression l\ ~ Z p £ 4 valid for single filaments, where in both 



cases l p = 17 pm |21[ denotes the persistence length of 



a single actin filament. Combining the above results and 
setting the deformation mode length A equal to l c one 
finally arrives at the following prediction 



Go ~ R z c w a ■ 6- 
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where the exponents are given by 2 = 2y — Ax and 
w = 7/3. Thus the scaling exponent of the plateau 
modulus can be related to parameters describing the mi- 
crostructure such as the scaling of the mesh size as well 
as the dependence of the bundle thickness and elastic- 
ity on R. From our measurements of a; = 0.27 and 
z = 1.5, and by assuming S to be a constant, a value of 
y = 1.29 is obtained, which is in reasonable agreement 
with the expected y « 1. This result is largely insensi- 
tive to the assumption of constant 5, since by assuming 
5 to change according to simple Langmuir kinetics an 
exponent y f=a 1.21 is obtained. 

To further characterize the elastic properties in the 
bundled regime, the non-linear elasticity of the network 
is investigated. For samples with R > R* a constant 
shear rate is applied and the resulting stress is reported. 
From the smoothed (7(7) relation the numerical deriva- 
tive yielding the differential modulus K = da/d-f is cal- 
culated (Fig. [3|). For small strains of 7 = 1 — 10 % 
linear behavior is observed, where the differential modu- 
lus follows K ~ R 15 in agreement with our oscillatory 
measurements. A non-linear response is observed above 
7c which is determined as the strain at which K de- 
viates by 5 % from its value in the linear regime. Up 
to R = 0.1 a strain hardening occurs while for very 
high values of R the linear regime is directly followed 
by strain weakening. The disappearance of the strain 
hardening at high concentrations of fascin might be the 
result of the rupturing of fascin-actin bonds - vcrysimilar 
to what was reported for rigor- HMM- networks Ml. The 
floppy mode description also has implications on the on- 
set of the non-linear behavior. As has been argued in [l(| 
large strains necessarily lead to stretching even if the de- 
formations were only of bending character. The bundle 
stretching A is related to the transverse displacement S n& 
by simple geometric considerations as Z 2 +<5 2 a = (l c +A) 2 . 
The floppy mode description thus applies as long as this 
stretching is small compared to the available thermal ex- 
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FIG. 3: Differential modulus K — da/d'y plotted versus de- 
formation 7 for fascin networks in the bundle phase (c Q = 
0.4 mg/ml and increasing R: diamonds R = 0.02, upright 
triangles: R — 0.05, circles: R = 0.1, downright triangles: 
R = 0.2, stars: R = 0.5). The inset shows the critical 
strain -y c in dependence on R. 



cess length A A ~ l c b/Nl p jl9|. This defines a critical 
strain 7c ~ l^b/N) 1 / 2 ~ R-y+x b -l/6 ~ iT 1 ' ^ 1 / 12 for 
the onset of non-linear effects. The scaling with R is in 
excellent agreement with our measurements (see inset of 
Fig. [3]). The weak dependence on the cross- linker spac- 
ing 7 C ~ <5 -1 / 12 also implies that this result is insensitive 
towards any putative dependence on the fascin concen- 
tration via S = S(cf). 

On the other hand, if one were to apply an affinc 
stretching model, a different picture emerges, where S(cf) 
has to be tuned to obtain a reasonable data fit. In such 
a model one would assume the modulus to be given by 
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where fc|| ~ a 3/2 (NKf) 2 jl\ is the stretch- 
ing stiffness of the bundle [ljj ]. The deformations arc 
assumed to be affine, implying (5 a ff ~ il c - The modulus 
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Supplement! ! while the 



thus reads as G a g ~ R 2x c a ■ S~ 

critical strain 7 C ~ 5 x / 2 R~ 2x is obtained by equating AA 
with the affinc deformation (5 a ff ~ "fl c . This model can 



-0 



with an exponent [3 



only fit the data by assuming S 

in the range of /3 = 0.6 — 0.9, such that for f3 = 0.6 
the i?-dependence of the modulus and for (3 = 0.9 the Co- 
dependence of the modulus is reproduced. To finally de- 
cide whether or not the application of an affine stretching 
model is equally successful as the floppy mode approach, 
5{cf) would have to be determined by scattering experi- 
ments. 

In summary, on the basis of a combined microscopy 
and rheology study we have shown that the actin bind- 
ing protein fascin mediates a transition from an entan- 
gled polymer solution to a homogeneously cross-linked 
bundle phase. These phases differ both in structure and 
mechanical properties. The location of the transition is 
given by a simple relation between the ABP concentra- 



tion and the entanglement length. Moreover this transi- 
tion point seems to be more general since it also occurs 
at similar ABP concentrations in other systems such as 
isotropically cross-linked networks or even composite net- 
works 0, uL [ItJ . The transition is a consequence of the 
interplay between the chemical kinetics of the binding 
proteins, the bending rigidity of the polymers and the 
entropic forces between those components. How the con- 
certed action of those driving forces leads to such a struc- 
tural transition is an interesting theoretical problem. The 
elasticity in the bundled phase is well explained in terms 
of a recently developed floppy mode picture [Io| . We 
have argued that in the absence of a significant amount 
of stored length in the bundles, non- affine bending is 
the dominant low energy excitation. It explains both 
the linear elasticity and the onset of non-linear behav- 
ior. This model has to be seen as an alternative to 
affine models where the elastic response is due to pulling 
out stored length fluctuations. While the elasticity of 
isotropic networks may be predominantly determined by 
such an entropic stretching of single filaments between 
the cross-linking points, we suppose that the elastic re- 
sponse of composite phases may instead be dominated 
by non-affinc deformations of bundles as described by 
floppy modes. The detailed understanding of the pre- 
sented purely bundled network, composed of shear dom- 
inated bundles, provides a benchmark for addressing the 
further challenge to describe the mechanics of networks, 
which are dominated by structural heterogeneities. 
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